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I.  SOFTWARE  PROCESSORS 


1.1  Introduction 


This  report  discusses  improved  software  bionic  sonar  imple- 
mentations (Volume  1),  spectrogram  analysis  (Volume  2),  and  the 
development  of  some  new  concepts  that  have  been  suggested  by  a 
study  of  animal  sonar  systems  (Volume  3). 

The  software  implementations  in  Volume  1 have  followed  two 
different  philosophies  for  sonar  target  recognition.  The  first  philosophy 
is  based  upon  the  idea  that  the  various  components  in  an  echo  can  be 
theoretically  predicted  from  the  geometrical  theory  of  diffraction, 
and  that  echoes  can  be  classified  by  detecting  these  components  and 
by  determining  their  relative  times  of  occurrence.  This  philosophy 
has  led  to  the  so-called  "bionic  component"  method. 

The  second  philosophy  for  software  target  classification  is 
to  empirically  determine  the  important  target  components  by  ana- 
lyzing the  echoes  from  a variety  of  objects.  The  important  com- 
ponents can  be  ascertained  by  means  of  principal  component 
analysis. 

The  biological  aspect  of  the  "bionic  component"  method 
arises  because  the  various  theoretical  components  are  best  detected 
with  a bank  of  proportional  bandwidth  filters  when  a dolphin-like 
pulse  is  used  as  a signal.  There  is  evidence  that  some  mammalian 
auditory  processes  can  be  modelled  with  such  a filter  bank. 

The  biological  aspect  of  the  principal  component  method 
arises  from  the  domain  in  which  echoes  are  analyzed  and  components 
are  derived.  The  peripheral  auditory  system  apparently  forms  a 
spectrogram-like  representation  of  a sonar  echo.  Another  domain 
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for  echo  representation  is  suggested  by  the  relevance  of  a Fourier- 
Mellin  transform  to  models  of  hearing.  Principal  component  classifiers 
for  both  spectrogram  and  Fourier -Mellin  echo  representations  have 
been  implemented  in  software. 

Three  different  software  algorithms  for  bionic  echo  classi- 
fication (bionic  component  detection,  principal  component  spectro- 
gram analysis,  and  principal  component  Fourier-Mellin  transform 
analysis)  have  been  tested  against  one  another  by  computing  mis- 
class  i '-ation  probabilities  in  noise  and  bottom  reverberation. 

These  tests  indicate  that  the  spectrogram  representation  gives  the 
best  performance,  with  a probability  of  misclassification  < 0.  1 
occurring  for  SNR>  0 dB  and  SCR  (signal-to-clutter  ratio)  5-  6 dB. 

The  superior  experimental  performance  of  a spectrogram 
correlator  indicates  tha*  this  method  should  be  analyzed  further 
from  both  optimum  detector  and  biological  modelling  viewpoints. 

This  analysis  is  the  subject  of  Volume  2. 

It  is  demonstrated  in  Volume  2 that  spectrogram  correlation 
is  a locally  optimum  detection  operation  for  low  SNR  signals  that 
have  been  passed  through  a random,  time-varying  filter  (e.  g.  , a 
range-extended  sonar  target  and/or  a long-range  or  shallow-water 
propagation  channel),  in  a Gaussian  noise  background.  For  non- 
Gaussian  noise,  the  local  optimality  of  the  process  is  preserved  by 
changing  the  power  law  that  is  associated  with  the  envelope  detection 
operation,  i.  e.  , by  correlating  stored  data  with  a modified 
spectrogram. 

Detection  and  estimation  performance  of  the  spectrogram 
correlation  process  is  predicted  and  compared  with  behavioral  data 
in  Volume  2.  Green's  summation  law  for  human  detection  of  a 
sum  of  gated  sinusoids  is  consistent  with  the  expected  performance 
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of  a spectrogram  correlator,  and  Siebert's  description  of  the 
behavioral  standard  deviation  of  an  auditory  frequency  estimate  in 
humans  is  also  consistent  with  a spectrogram  correlator  model. 
Neither  of  the  above  behavioral  results  is  consistent  with  a matched 
filter  model. 

In  the  "new  concepts"  category,  the  following  results  have 
been  obtained,  and  are  discussed  in  detail  in  Volume  3: 

1.  A new,  range-angle  ambiguity  function  has  been  formu- 
lated in  order  to  combine  signal  and  array  optimization  for  high  reso- 
lution sonar.  The  new  ambiguity  function  helps  to  demonstrate  how 
increased  system  bandwidth  can  offset  structural  array  constraints 
(two  ears,  separated  by  a small  distance)  in  animal  sonar  systems. 

2.  A feasibility  study  has  indicated  that  bats  and  dolphins 
may  be  able  to  estimate  cross -range  velocity  when  targets  are  less 
than  1 m.  away.  This  capability  may  provide  the  animals  with  a 
counter-countermeasure  to  cope  with  evasive  maneuvers.  In  a man- 
made system,  the  measurement  of  cross-range  velocity  at  longer 
distances  can  be  accomplished  with  more  powerful  sonars  that  have 
larger  time -bandwidth  products,  coherent  pulse-to- pulse  integra- 
tion, or  an  array  that  is  larger  than  an  animal's  interaural  distance. 

3.  A test  is  proposed  to  distinguish  a locally  optimum 
version  of  the  energy  spectrum  analyzer  of  R.  A.  Johnson  from  a 
matched  filter  model  for  animal  echolocation. 

In  the  remainder  of  this  volume,  three  different  software 
processors  or  computerized  models  for  animal  sonar  systems  are 
discussed  and  compared.  The  discussions  include  a justification  of 
each  model  in  terms  of  biological  results,  as  well  as  a description 
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of  the  signal  processing  operations  that  are  performed.  Target 
classification  performance  of  the  models  is  then  compared  for 
various  levels  of  background  noise  and  bottom  reverberation. 

1.2  The  "Bionic  Component"  Model 

The  stimulus  features  that  are  important  to  an  animal  can 
be  hypothesized  by  observing  the  differences  in  signals  that 
would  be  useful  for  stimulus  classification.  In  the  case  of 
echolocation,  the  stimulus  is  an  echo  and  stimulus  classification 
depends  upon  differences  between  echoes  from  various  targets. 

By  using  modern  pattern  recognition  techniques,1  it  is  possible  to 
empirically  determine  important  classification  features  by  obtain- 
ing a large  data  base  of  echoes  from  objects  in  the  animal's  environ- 
ment. It  is  also  possible  to  use  theoretical  predictions  of  target 
scattering  characteristics  in  order  to  decide  which  echo  features 
should  be  most  important.  The  "bionic  component"  technique  uses 
the  latter  philosophy. 

Research  on  the  prediction  of  target  scattering  characteris- 
tics has  resulted  in  two  main  theoretical  approaches.  Geometrical 

and  physical  optics  have  evolved  into  a technique  that  is  known  as 

0.0 

the  geometrical  theory  of  diffraction.  This  theory  involves  the 
prediction  of  scattering  characteristics  from  separate  structures 
such  as  the  tip  of  a cone  or  an  edge,  and  it  includes  creeping  wave 
effects  (i.e.,  waves  that  propagate  around  the  surface  of  the 
object).  The  echo  components  from  various  target  structures  are 
combined  with  proper  delay,  in  order  to  form  the  composite  echo. 
The  geometrical  theory  of  diffraction  is  primarily  a high  frequency 
approach  for  wavelengths  that  are  smaller  than  the  target. 
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The  second  approach  to  target  scattering  prediction  is  the 
use  of  natural  modes  or  resonances,  as  derived,  for  example,  with 
the  singularity  expansion  technique.  Natural  modes  tend  to  domi- 
nate the  trailing  edge  of  an  echo,  whereas  specular  components 
(mirror-like  reflections  or  glints)  with  higher  amplitudes  are 
found  at  the  beginning  of  the  echo.  The  advantage  of  natural  mode 
characterization  is  that  the  resonant  frequencies  are  independent 
of  aspect,  and  it  would  thus  seem  that  the  target  echo  need  not  be 
modelled  as  an  aspect-dependent  random  process.  Experimental 
results  with  various  sonar  targets,  however,  have  shown  that, 
although  natural  resonances  may  be  independent  of  aspect,  the  exci- 
tation of  a given  resonance  and  the  energy  that  is  re-radiated  toward 
the  receiver  is,  in  fact,  aspect-dependent.  Furthermore,  the 
trailing  edge  of  the  echo  tends  to  be  attenuated  relative  to  the 
specular  parts  of  the  echo.  Finally,  the  trailing  edge  of  an  echo 
is  difficult  to  analyze  if  the  target  is  immersed  in  clutter. 

According  to  Bechtel  and  Ross,8  the  geometrical  theory  of 

diffraction  predicts  that  many  scattering  centers  will  have  cross- 

sections  that  vary  as  Xn,  where  X is  wavelength  and  -2  ?n  ? 2, 

See  Table  1.  If  cross  section  varies  as  An,  then  the  associated 

range-dependent  transfer  function  must  vary  as  a half-power  of 

n / 2 

frequency,  and  the  echo  spectrum  is  (jw)  U(w)  for  -2  < n < 2, 

where  U(cj)  is  the  Fourier  transform  of  the  transmitted  signal.  An 

idealized  highlight  echo  can  be  synthesized  by  passing  an  impulse 

n/2 

through  a filter  with  transfer  function  (jcu)  U(u>).  The  parameter 
n is  range  dependent,  and  changes  from  one  highlight  to  the  next. 

(A  highlight  is  a local  maximum  in  the  echo  envelope.  ) Since  tar- 
gets are  distributed  in  azimuth  and  elevation  as  well  as  in  ra  >g-», 
it  is  possible  that  more  than  one  type  of  geometry  in  Table  1 will 
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Table  1.  Geometrical  diffraction  theory  applied  to  scattering  centers 
(monostatic  case).  * 


From  Bechtel  and  Ross,  Cornell  Aeronautical  Lab.  Report  ER/RIS-10. 


be  encountered  in  a given  range  interval.  In  this  case,  an  idealized 
target  is  modelled  by  a sum  of  two  or  more  different  transfer 
functions. 

Recent  experimental  results14  seem  to  indicate  that  a high- 
light is  rarely  comprised  of  an  echo  from  a single  point  in  range. 
Sonar  highlights  are  caused  by  range-distributed  impedance  changes. 
An  echo  from  a range-distributed  structure  can  be  simulated  by 
passing  time- gated  shot  noise  (closely  spaced  delta  functions  with 
uncorrelated  epochs)  through  a filter  that  has  a range-dependent 
transfer  function.  In  terms  of  the  geometrical  theory  of  diffraction, 
highlight  echoes  are  best  modelled  by  passing  time -gated  shot  noise, 
rather  than  a single  impulse,  through  a filter  with  transfer  function 
(jo))n^2U(a >),  where  -2  < n < 2 and  U(w)  is  the  spectrum  of  the  trans- 
mitted signal.  This  random  aspect  of  the  model  leads  to  a receiver 
that  is  different  from  the  original  "bionic  sonar"  proposed  by 
Altes.15 

A method  of  target  characterization  that  is  based  upon  the 

above  model  is  to  implement  a detector  that  decides  whether  a given 

echo  component  is  present  at  any  given  range.  Since  each  echo 

component  is  modeled  as  a sample  function  of  a random  signal  with 

known  power  spectral  density,  the  detection  can  be  implemented 

16 

with  an  estimator -correlator  process. 

The  estimation  part  of  the  process  consists  of  passing  echo 
data  through  each  of  five  Wiener  filters.  Since  any  given  echo 
component  constitutes  an  independent,  interfering  signal  for  the 
reception  of  a different  component,  the  filters  tend  to  be  narrow- 
band,  and  the  m^  filter  transfer  function  is  centered  near  the 
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maximum  point  of  the  power  spectral  density  |(jaj)  U(u>)|  of  the 
m**1  echo  component.  Idealized  versions  of  the  five  filters  can  be 
obtained  by  slicing  the  overlapping  echo  components  at  their  cross- 
over levela  in  the  frequency  domain,  as  shown  in  Figure  1.  The 
idealized  estimating  filters  are 

Wm (w)  = max  [|(jt*>)m/2U(<J)|2  - c,  o]  , 

m = -2,  -1,  0,  1,  2 (1) 

where  c is  the  cross-over  level  that  is  shown  in  Figure  1. 

The  correlation  part  of  the  process  is  accomplished  by 
multiplying  each  filter  output  by  the  original  data  and  integrating 
over  the  expected  duration  of  the  highlight  echo.  The  integrator 
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Figure  1.  Idealized  Wiener  filters  Wm(o>)  are  defined  as 
= max{ii)m|U(u>|2  “ c,  OJ  where  c is 
the  cross-over  level  fV»r  the  amplitude  normal- 
ized echo  components  with  power  spectral 
densities  wm|U(<j)|2,  m = -2,  ....  2. 
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output  must  exceed  two  separate  thresholds  in  order  to  trigger  a 
decision  that  the  component  is  present  at  any  given  range. 

The  first  threshold  is  determined  by  expected  echo  energy. 
The  second  threshold  varies  with  range,  and  is  set  to  nine-tenths 
of  the  largest  estimator-correlator  response  at  any  given  time. 

The  purpose  of  the  second  threshold  is  to  compensate  for  the 
cross-talk  or  excitation  of  the  m^  filter  by  the  n**1  echo  compo- 
nent. From  Figure  1,  it  is  apparent  that  the  n^  component  can 
excite  the  m^  filter,  even  if  the  filter  transfer  functions  do  not 
overlap. 

Application  of  the  detection  algorithm  to  a sonar  echo 
results  in  a pattern  of  suprathreshold  responses,  as  shown  in 
Figure  2.  Different  patterns  are  obtained  for  different  targets, 
and  the  patterns  resemble  a one-bit  quantization  of  a spectrogram. 

If  the  signal  has  log-normal  spectral  amplitude,  the  filters 
W (w)  in  (1)  have  bandwidths  that  are  proportional  to  frequency. 

The  existence  of  proportional  bandwidth  filters  in  mammalian  audi- 
tory systems  can  be  inferred  from  the  results  of  noise-masking 
experiments  and  from  neurophysiological  data.  Many  ceta- 

cean echolocation  signals  have  spectral  magnitudes  that  are  approxi- 
mately log-normal  in  shape.15 

In  a previous  report,14  only  the  total  number  of  supra- 
threshold responses  from  each  filter  was  used  for  target  classifi- 
cation. This  method  disregards  the  relative  timing  information 
that  appears  in  Figure  2.  It  was  therefore  decided  to  add  a timing 
parameter  that  describes  the  mean  time  for  all  the  suprathreshold 
responses  relative  to  the  mean  time  for  the  filter  that  has  the  largest 
number  of  responses.  This  timing  parameter,  along  with  the  num- 
ber of  detections  from  each  filter,  resulted  in  six  target  features 
altogether.  The  timing  parameter  is  computed  as  follows. 
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Let  z..  be  the  response  of  the  j**1  detector  at  time  t..  If 
lj  1 

the  response  exceeds  threshold,  then  z..  - 1,  and  z..  = 0 otherwise. 

lJ  1J  fh 

The  mean  time  for  the  suprathreshold  responses  from  the  j filter 
is  then 


T. 

J 


t.  z. . 
* U 


/ E 


Z..  . 

IJ 


Suppose  that  the  maximum  number  of  detections  is  observed  at  the 

output  of  detector  number  J.  The  mean  time  of  all  the  suprathreshold 

detector  responses  relative  to  T is  then 

J 

5 5 


(2) 


The  timing  parameter  AT  in  (2)  is  an  extra  feature  that 
was  used  to  generate  new  misclassification  probabilities  for  the 
bionic  component  method.  The  performance  of  this  algorithm  will 
be  compared  with  that  of  the  other  two  models  in  Section  1,  5.  The 
five- component  system  without  AT  will  also  be  included  in  the 
comparison. 
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1.3 


The  Spectrogram  Model 


A spectrogram  is  a display  of  a sequence  of  short  duration 

energy  density  spectra,  and  it  portrays  the  distribution  of  signal 

energy  in  both  frequency  and  time.  A spectrogram-like  signal  repre- 

21-23 

sentation  may  be  formed  by  the  peripheral  auditory  system.  The 

use  of  the  spectrogram  in  sonar  echo  processing  depends  upon  the 
integration  time,  i.  e. , the  duration  of  the  time  window  or  the 
impulse  response  duration  of  the  filters  that  are  used  to  generate 
a sequence  of  energy  density  spectra.  When  the  spectrogram  is 
formed  by  a bank  of  filters  with  envelope  detected  responses,  the 
integration  time  varies  inversely  with  the  bandwidth  of  each  filter. 
Unfortunately,  there  is  a large  spread  of  possible  bandwidths  and 
integration  times  that  have  been  measured  by  noise-masking  experi- 
ments, frequency  discrimination  studies,  and  neurophysiological 
investigations.  It  would  appear,  in  fact,  that  bandwidth  (and  inte- 
gration time)  depends  upon  the  particular  detection  problem  that 

24  25 

is  presented  to  a subject.  ’ 

If  long  integration  times  are  relevant,  then  both  a trans- 
mitted pulse  and  the  resulting  echoes  are  included  in  the  Fourier 

26 

transform,  and  the  Johnson- Titlebaum  model  seems  to  be  appro- 
priate. For  very  short  integration  times,  one  can  divide  up  a target 
into  a sequence  of  range  cells,  and  observe  the  evolution  of  reflected 
spectral  energy  as  a function  of  range.  Both  long  and  short  integra- 
tion times  have  certain  advantages,  but  we  consider  only  a short 
integration  time  in  this  volume. 

The  window  length  for  our  software  algorithm  is  9 cm  long 
in  range,  and  a new  window  is  generated  every  1.8  cm.  One  aspect 
of  the  software  spectrogram  processor  may  not  be  particularly 
"bionic."  Before  echoes  are  time -windowed  and  Fourier  transformed. 
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they  are  passed  through  a filter  that  is  matched  to  the  transmitted 
signal.  Since  the  time -bandwidth  product  of  the  signal  (a  typical 
dolphin  pulse)  is  not  very  large,  this  operation  is  probably  optional. 
For  test  purposes,  however,  it  was  decided  that  matched  filtering 
would  improve  the  temporal  resolution  of  the  spectrogram. 

A Fortran  program  for  computing  a spectrogram  for  wide- 
band sonar  echo  analysis  is  listed  in  Appendix  A. 

It  has  been  mentioned  that  an  alternative  to  theoretical 
scattering  prediction  is  to  use  empirical  echo  data  in  order  to 
obtain  a set  of  features  for  target  classification.  The  weakness  of 
such  an  approach  is  that  the  data  base  that  is  used  for  feature  defi- 
nition may  be  too  limited.  On  the  other  hand,  the  use  of  wideband 
signals  upon  targets  that  are  made  of  different  materials  can 
severely  degrade  the  reliability  of  theoretical  predictors.  Short- 
term spectrogram  analysis  has  therefore  been  coupled  with  empirical 
feature  determination. 

Features  are  significant  attributes  of  the  data  that  can  be 
used  for  classification.  Consider  a filter  that  processes  the  data 
in  order  to  extract  a particular  feature.  A set  of  such  filters  with 
orthogonal  impulse  responses  would  constitute  an  orthogonal  basis 
for  decomposition  of  the  data.  Since  the  desired  number  of  features 
is  limited  in  practice,  the  orthogonal  basis  set  is  generally  incom- 
plete, i.  e. , it  is  not  capable  of  reconstructing  every  possible  data 
waveform.  One  can  obtain  a "best"  incomplete  N-dimensional  basis 
set  by  finding  the  set  that  gives  the  least  expected  mean- square  error 
between  the  reconstructed  data  and  a representative  set  of  the  ori- 
ginal data,  for  a given  number  of  feature  vectors  (orthogonal  filters), 

N.  The  resulting  set  of  N feature  vectors  are  the  N eigenvectors 

27 

of  the  data  covariance  matrix  that  have  the  largest  eigenvalues. 
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The  associated  eigenvalues  are  indicative  of  the  amount  of  informa- 
tion that  is  contributed  by  each  feature  vector,  i.  e. , the  percentage 
of  the  total  data  variance  that  is  accounted  for  by  each  component. 

Application  of  the  "best"  incomplete  basis  set  corresponds 
to  a Karhounen- Loeve  decomposition  of  the  data,  and  it  is  also 
known  as  a principal  component  analysis.1  Since  a limited  number 
of  basis  functions  are  used  to  convey  the  maximum  possible  infor- 
mation about  the  data,  one  can  assume  (in  the  absence  of  a better 
approach)  that  these  basis  functions  correspond  to  the  most  important 
features  for  classification 

Typical  spectrogram  responses  to  various  echoes  are  shown 

in  Figure  3.  A sampled  two-dimensional  plot  can  be  represented 

as  a long  column  vector  by  sequentially  listing  the  elements  in  the 

first  row,  then  the  second  row,  etc.  A column  vector  of  this  kind 

was  constructed  from  the  echo  spectrogram  of  each  of  the  37  targets 

in  Table  2.  The  resulting  set  of  column  vectors  was  then  used  to 

construct  a data  covariance  matrix,  C . If  d . is  the  column  vector 

s — 1 

that  is  derived  from  the  spectrogram  of  the  i**1  echo,  then 

37 

c„  - if  Z)  (4i  - 2)  <3,  - I>T  <3> 

i=  1 

where  the  T denotes  transpose  and 


37 


i=  1 


The  target-induced  variation  in  the  magnitude  of  the  spectrogram 
at  each  point,  i.  e. , the  variation  that  is  caused  by  different  echoes, 
thus  contributes  to  the  derivation  of  a feature  set. 
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Figure  3b.  Sampled  spectrogram  of  the  match-filtered  echo  from  a 
solid  aluminum  sphere,  3"  diameter. 
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Figure  3c.  Sampled  spectrogram  of  the  match-filtered  echo  from  a 
hollow  aluminum  cylinder,  1. 5"  diameter  x 7"  long  x 
3/l6"  thick,  broadside  aspect. 


Table  2.  List  of  targets  that  were  used  for  the  three 
software  processors. 


Unless  otherwise  noted,  targets  are  broadside  to  the  direction 
of  propagation.  Broadside  - 0° , end-on  = 90".  Cylinder 
dimensions  are  diameter  x length  x wall  thickness.  Hollow 
cylinders  are  open-ended  (water-filled). 


5"  biconic 
3"  sphere 
2.  5 x 7"  rock  cyl. 

1.5  x 7 x 3/16"  Al.  cyl. 

1.5  x 7 x 4/16"  Al.  cyl. 

1.5  x 7 x 5/16"  Al.  cyl. 

1.5  x 7 x 6/16"  Al.  cyl. 

3 x 7 x 2/16"  Al.  cyl. 

3x7  *4/16"  Al.  cyl. 

3 x 7 x 6/16"  Al.  cyl. 

3 x 7 x 8/16"  Al.  cyl. 

3 x 7 x 10/16"  Al.  cyl. 

3 x 7 x 8/16"  PVC  cyl. 

Transmitted  signal  (planar 

target) 

6 cm  foam  cube,  0° 

6 cm  foam  cube,  45° 

6 cm  x 6 cm  foam  cyl. , 0° 
6 cm  x 6 cm  foam  cyl.  , 90° 

2.5  x 7"  Al.  cyl.,  0° 
2.5x7  Al.  cyl. , 5° 

2.5  x 7"A1.  cyl.,  10° 

2.5  x 1 Al.  cyl.,  15* 


2.  5 x 7' 
2.  5 x 7' 
2.  5 x 7' 
2.  5 x 7' 
2.5x7' 
2.  5 x 7' 
2.  5 x 7' 
2.5x7' 
2.  5 x T 
2.5x7' 
2.5x7 
2.5x7 
2.  5 x 7 
2.  5 x 7 
2.  5 x 7 


Al.  cyl.,  20* 


Al.  cyl. 
Al.  cyl. 
Al.  cyl. 


, 25* 
, 30* 
. 35* 


Al.  cyl. , 40* 


Al.  cyl. 
Al.  cyl. 
Al.  cyl. 


. 55* 


Al.  cyl.,  60* 
Al.  cyl.,  65° 
Al.  cyl.  , 70* 


Al.  cyl. 
Al.  cyl. 
Al.  cyl. 


, 80* 
, 85° 


Al.  cyl.,  90° 


Table  3 (middle  column)  shows  the  eigenvalue  estimates 

for  the  first  sixteen  eigenvectors  of  C . 

s 

The  first  seven  eigenvectors  of  C are  shown  in  Figure  4. 

s 

These  functions  of  frequency  and  time  represent  the  features  that 
are  most  useful  for  the  classification  of  the  targets  in  Table  2 
from  a spectrogram  representation  of  wideband  sonar  echoes. 

The  functions  in  Figure  4 have  interesting  implications  for 
models  of  the  acoustic  and  visual  systems.  Suppose  that  each  of 
the  features  in  Figure  4 corresponds  to  a neuron  in  the  auditory 
cortex.  The  neuron  that  is  associated  with  Figure  4b,  for  example, 
has  an  inhibitory  center  at  40  kHz,  surrounded  by  an  excitatory 
area  in  frequency  and  time.  The  neuron  in  Figure  4c  is  excited 
by  spectrogram  outputs  that  appear  early  in  time,  and  it  is  inhibited 
by  outputs  that  appear  later  in  time.  Such  neurons  display  not 
only  lateral  inhibition  for  neural  "sharpening"  in  the  frequency 
domain,28  but  also  forward  and/or  backward  masking  in  the  time 
domain.29  The  similarity  to  receptive  fields  in  vision-30  is 

obvious.  Our  results  suggest  that  visual  cells  with  specialized 
receptive  fields  may  form  an  orthonormal  basis  for  pattern  recog- 
nition, and  the  well-known  excitatory  center /inhibitory  surround 
cell  is  one  of  the  principal  components  for  such  a basis. 

Before  principal  component  analysis  occurs,  a data  spectro- 
gram must  be  registered  in  time,  i.e.,  the  temporal  position  of  the 
data  relative  to  the  basis  set  must  be  determined.  The  basis  set 
covers  a frequency-time  window  that  has  dimensions  56  kHz  x 672 
psec.  If  the  match-filtered  echo  data  is  denoted  r(t.),  where  the 


Table  3. 

Sixteen 

C and 
s 

largest  normalized  eigenvalues 
Cfm  in  descending  order. 

of 

Eigenvalue  Number 

C 

^fm 

1 

.459 

. 284 

2 

. 202 

. 198 

3 

. 100 

. 13! 

4 

. 068 

. 095 

5 

. 042 

. 057 

6 

. 036 

. 044 

7 

. 031 

. 031 

8 

. 021 

. 025 

9 

. Oil 

. 022 

10 

. 009 

. 018 

11 

. 007 

. 015 

12 

. 005 

. 014 

13 

. 004 

. 012 

14 

. 003 

. 010 

15 

. 002 

. 008 

16 

. 001 

. 007 

zo 
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BOTTOM  VIEW 


Figure  4a. 

Figure  4.  Principal  components  or  basis  functions  for 

the  analysis  of  echo  spectrograms.  Figure  4a 
shows  the  component  with  the  largest  eigen- 
value; Figure  4g  shows  the  component  with  the 
seventh  largest  eigenvalue. 


sampling  interval  is  24  psec  long,  then  the  starting  time  for  spectro- 
gram analysis  (i.e.,  the  left-hand  edge  of  the  frequency-time  win- 
dow) corresponds  to  the  point  t.  where 

28 

sum  (t.)  = [r(tj+i-l)  ] 2 (5) 

i=l 

is  maximized.  The  above  sum  is  a local  energy  detector  fo*-  esti- 
mating the  range  of  an  extended  target  with  a &?2  p sec  long  echo, 
where  the  energy  is  computed  at  the  output  of  a filter  that  is  matched 
to  the  transmitted  signal. 

The  target  classification  performance  of  the  spectrogram  proc 
essor  has  been  determined  in  backgrounds  of  white  noise  and  bottom 
reverberation.  This  performance  is  compared  to  the  other  two 
software  processors  in  Section  1.  5. 

1,  4 The  Fourier-Mellin  Transform  Model 

The  Fourier-Mellin  transform  of  an  analytic  signal  u(t)  is 
defined  as 

Cu0*>  = f u(f,  «jxl°8(2,'f)df/f  (6, 

0 

where  U(f)  is  the  Fourier  transform  of  u(t).  This  transform  is 

relevant  to  biological  systems  because  the  phase  function  exp[  jx  log 

(2?rf)]  can  be  synthesized  with  a bank  of  proportional  bandwidth 

filters,31, 33  and  because  the  transform  can  be  used  to  obtain  a scale- 

and  shift- invariant  waveform  representation.  The  process  in  (6) 

is  equivalent  to  matched  filtering  of  a linear  period  modulated 

rectangular  pulse,  ’ i.  e.  , a waveform  that  is  used  by  some  bats 

33*35 

and  dolphins  for  echolocation. 
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A time-dependent,  two-dimensional  signal  representation 
can  be  obtained  from  (6)  by  defining  a set  of  band  limited  filters 


V«  * 


f exp[  log(2fff)]  , f <f<f 


0 , I < fL  or  l > fH 


where  k = 1,  2,  . . . , M. 

The  two-dimensional  representation  is  then 


rH 

Fuv(ti'V  = / U(f)  Vk(f)  ej27rfti  df 


A delayed  signal  u(t-r)  with  Fourier  transform 


U(f)  = 


f exp[  -j(2JTfr  + c log  (27Tf)]  , f < f < f 

L T1 


0 , f < fL  or  f >fH 


will  be  match -filtered  by  the  process  in  (8),  and  the  envelope- 

detected  matched  filter  response  will  be  observed  on  the  x = c 

k 

profile  of  the  function  F (t.,  x ).  The  maximum  matched  filter 

UV  X K 

response  will  occur  at  t.  = r,  x^  = c. 

Another  bionic  aspect  of  the  transformation  in  (8)  is  that 
the  filters  Vk(f)  are  matched  to  rectangular  pulses  with  different 

chirp  rates,  x^.  The  existence  in  humans  of  filters  that  are  chirp- 

or 

rate- sensitive  has  been  deduced  by  Kay  and  Mathews. 


The  function 


F„.<'-V  =1  / U(f)Vk(f)ej2"£,df|2. 


k = 1,  2,  . . . , M (10) 


where 


vk(f)  = 


f exp  [ j log  (2  7r  f ) ] , f C 


0,  otherwise 


will  be  called  the  FM-gram,  since  it  is  a modified  spectrogram 
that  is  based  upon  the  Fourier-Mellin  transform  rather  than  upon 
an  ordinary  spectral  representation. 

A Fortran  program  for  computing  the  FM-gram  is  listed  in 
Appendix  B.  Seven  k-values  were  used,  corresponding  to 


x = (k-3)  (-0.  25)  (2r)  / (log  1.2),  k=l,2,  ....  7.  (12) 


Examples  of  FM-gram  echo  representations  for  three  tar- 
gets are  shown  in  Figure  5.  Large  differences  between  the  three 
echo  representations  are  immediately  noticeable.  If  such  large 
differences  are  found  for  each  target,  then  one  would  expect  that 
the  various  target  representations  are  nearly  orthogonal,  and  that 
a large  number  of  principal  components  are  needed  in  order  to 
approximate  all  the  echoes  in  the  data  set.  This  situation  will  be 
indicated  by  a relatively  slow  drop-off  of  the  eigenvalues  of  the 
target  covariance  matrix,  C This  phenomenon  can  indeed  be 
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Figure  5a.  FM-gratn  samples  for  the  echo  from  a solid 
aluminum  biconic  (frustrum  of  a cone) . Largest  response 
(98)  occurs  at  output  of  Filter  No.  4 (see  next  page). 
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Figure  5a.  (Continued) 
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Figure  5b.  FM-gram  samples  for  the  echo  from  a solid  aluminum 
3"  sphere.  Largest  response  (98)  occurs  at  output  of 
Filter  No.  7. 
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Figure  5c.  FM-gram  samples  for  the  echo  from  a hollow 
aluminum  cylinder  (1.5"  x 7"  x 3/16"  wall  thickness). 
Largest  response  (98)  occurs  at  output  of  Filter  No.  3. 
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observed  in  Table  3.  The  biological  counterparts  of  the  eigen- 
vectors of  C_.  . would  be  neural  receptive  fields  that  are  inhibited 
FM 

or  excited  by  frequency  modulated  pulses  with  various  chi.p  rates 
and  relative  times  of  occurrence. 

The  target  classification  performance  of  the  FM-gram  proces- 
sor has  been  determined  in  backgrounds  of  white  noise  and 
bottom  reverberation.  This  performance  is  compared  to  the  other 
software  processors  in  the  next  section. 

1.5  Comparative  Performance  of  Software  Processors 

The  performance  of  the  three  bionic  target  classifiers  was 
determined  by  grouping  targets  into  various  classes  and  by  com- 
puting probability  of  misclassification  in  a background  of  clutter 
(bottom  reverberation)  or  white  noise.  Probability  of  misclassifi- 
cation  for  various  signal-to-noise  ratios  (SNR)  and  signal-to- 
clutter  ratios  (SCR)  were  determined  by  adding  clutter  or  noise  to 
the  target  echoes  and  then  processing  the  resulting  data  with  each 
of  the  three  software  algorithms.  SNR  and  SCR  are  defined  as 
follows. 

SNR  = echo  energy/  average  noise  power  (13) 

where  the  noise  power  is  measured  at  the  output  of  a bandpass 
filter  with  low  frequency  cutoff  f = 5 kHz  and  high  frequency  cut- 
off f =95  kHz.  This  filter  simulates  the  pass  band  of  the  experi- 
H 

mental  sonar  system  that  was  used  to  gather  the  data.  Similarly, 

SCR  = echo  energy  / average  clutter  power  (14) 


— --- 


35 


where  the  clutter  is  a sample  function  of  bottom  reverberation.  The 
reverberation  was  measured  with  a sonar  signal  that  had  a higher 
chirp  rate  than  the  signal  u(t)  that  was  used  in  obtaining  target 
echoes.  The  reverberation  data  was  therefore  pre-processed  in 
software  in  order  to  obtain  clutter  that  would  have  been  measured 
with  u(t)  as  the  transmitted  signal.  A large,  anomalous  low 
frequency  clutter  component  was  observed  at  2 kHz,  and  this  com- 
ponent was  filtered  out  by  passing  clutter  echoes  through  a band- 
pass filter  with  f = 5 kHz  and  f = 95  kHz.  The  clutter  power  is 
L H 

measured  at  the  output  of  this  bandpass  filter. 

The  eigenvalues  of  Cg  in  Table  3 decrease  sufficiently 
rapidly  that  adequate  performance  can  be  obtained  with  only  seven 
eigenfunctions.  A classification  test  with  seven  eigenfunctions  of 
C gives  comparatively  poor  results.  The  reason  for  this  rela- 
tively poor  performance  is  that  target  information  is  more  distri- 
buted for  the  FM-gram  representation,  and  more  basis  vectors  are 
needed  in  order  to  adequately  describe  each  target.  This  effect  is 
represented  by  the  relatively  even  distribution  of  eigenvalue  ampli- 
tudes for  C in  Table  3.  The  classification  performance  of  an 
FM 

FM-gram  representation  has  therefore  been  evaluated  by  using 
sixteen  eigenfunctions  rather  than  the  seven  that  are  used  for  the 
spectrogram  representation. 

Is  it  advantageous  to  "spread  out"  the  target  information 
among  a comparatively  large  number  of  orthogonal  components,  or 
should  one  always  use  the  representation  that  has  the  least  number 
of  basis  vectors?  The  latter  strategy  corresponds  to  dimensionality 
reduction  in  a pattern  recognition  problem.  The  strategy  of  "spread- 
ing out"  information  corresponds  to  a high-diversity  system  as 
employed,  for  example,  in  spread- spectrum  communication.  By 
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using  the  Kullback- Leibler  information  number,  R.  F.  Daly37  has 
shown  that  an  optimum  degree  of  diversity  is  such  as  to  make  the 
signal-to-noise  ratio  for  each  orthogonal  component  equal  to  approxi- 
mately 2.  Daly's  result  suggests  that  better  performance  will  be 
obtained  at  low  SNR  if  the  available  echo  energy  can  be  concentrated 
among  a few  principal  components,  as  in  the  spectrogram  representa- 
tion. As  SNR  is  increased,  a more  diverse  echo  representation  is 
permissible.  We  might  therefore  predict  that  the  seven-component 
spectrogram  analysis  will  be  superior  to  FM-gram  with  sixteen 
principal  components  for  low  SNR.  It  will  be  seen  that  this  prediction 
is  commensurate  with  our  experimental  results. 

Target  echoes  were  classified  by  means  of  a mean- square 
distance  measure 


N 


where 


D. 

1 


r = 


(r  - f .) 
n ni 


(15) 


n=  1 


f . 
m 


N 


correlation  between  the  spectrogram  or  FM-gram  and 
the  n**1  eigenvector,  or  the  response  of  one  of  the  five 
suprathreshold  detection  counters  or  the  timing  param- 
eter for  the  "bionic  component"  method. 

n^1  eigenvector  correlator  or  "bionic  component" 
th 

response  for  the  i target  in  the  absence  of  noise  or 
reverberation  (archetypes) 

5 for  the  bionic  component  method  without  the  timing 
parameter 

6 for  the  bionic  component  method  with  the  timing 
parameter 
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N =7  for  principal  component  spectrogram  analysis 
= 16  for  principal  component  FM-gram  analysis 

If  D < D for  all  n i i,  then  the  echo  is  classified  as  target  i.  A 
i n' 

misclassification  occurs  when  the  test  (15)  judges  a target  to  be  a 
member  of  class  B when  the  target  has  been  assigned  to  class  A, 
or  vice-versa.  If  target  i is  wrongly  identified  as  target  j,  no 
misclassification  results  if  j is  in  the  same  class  as  i. 

Probabilities  of  misclassification  as  a function  of  input 
SNR  or  SCR  are  shown  in  Figures  6 and  7.  Figures  6a  and  7a  are 
concerned  with  distinguishing  a 2.  5 x 7"  solid  aluminum  cylinder, 
as  seen  from  various  aspects,  from  a variety  of  other  targets, 
which  are  listed  in  Table  2.  Figures  f>b  and  7b  display  misclassifi- 
cation probabilities  for  a cylinder  versus  non-cylinder  classifier, 
i.e.,  the  targets  in  Table  2 are  categorized  as  "cylinders"  or 
"something  else.  " Figures  6c  and  7c  show  results  for  a metal 
versus  non-metal  classifier,  again  using  the  targets  in  Table  2. 

The  graphs  in  Figures  6 and  7 include  probabilities  that 
would  be  obtained  by  random  guessing.  These  probabilities  are 
computed  from  the  equation 

Pr(misclass'n)  = Pr (misclass 'n| A)  Pr(A) 

+ Pr(misclass'n  |B)  Pr(B) 

where  the  target  set  is  divided  into  two  classes  A and  B.  Pr(A)  is 
the  probability  that  a target  from  class  A is  encountered  and 
Pr (misclas s 'n  |A)  is  the  probability  of  misclassification  given  that 
a target  from  class  A is  the  correct  choice,  and  a random  guess  is 
made.  For  example,  if  A is  the  class  of  32  cylinders  and  B is  the 
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Figure  6.  Relative  performance  of  software  algorithms  in 
additive  white  noise. 
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Figure  7.  Relative  performance  of  software  algorithms  in 
additive  reflections  from  bottom  clutter. 
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class  of  5 non-cylinders  in  Table  2,  then 

Pr(misclass’n)  = (^)(jf)  * (if)  (^)  = °' 23‘ 

The  following  conclusions  can  be  drawn  from  Figures  & 

and  7: 

1.  Spectrogram  correlation  using  seven  principal  components 
is  superior  to  the  other  algorithms  that  were  tested  for 
Pr  [misclassification]  < 0.  1.  For  spectrogram  analysis, 
probability  of  misclassification  drops  below  0.  1 for  SNR 
greater  than  -5  dB  to  0 dB  for  a white  noise  background,  and 
for  SCR  greater  than  1 dB  to  6 dB  for  a background  of  bottom 
clutter. 

2.  The  addition  of  a timing  parameter  to  the  bionic  component 
(or  geometrical  theory  of  diffraction)  method  is  helpful  for 
a white  noise  background  for  SNR  > 0 dB,  but  the  timing 
parameter  makes  little  difference  in  clutter  or  in  white 
noise  with  SNR  < 0 dB.  Performance  of  the  five -component 
GTD  method  seems  to  be  roughly  the  same  in  noise  and  rever- 
beration, but  the  spectrogram  method  does  worse  in 
reverberation  than  it  does  in  white  noise. 

3.  The  Fourier -Mellin  representation,  FM-gram,  has  larger 
storage  requirements  than  the  spectrogram  method,  since 
more  principal  components  must  be  used.  The  relative  per- 
formance of  FM-gram  in  white  noise  is  quite  poor,  but  the 
performance  is  better  in  clutter.  In  any  case,  FM-gram  is 
generally  inferior  to  the  spectrogram  method. 
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Since  the  spectrogram  correlator  with  seven  principal  com- 
ponents was  found  to  be  superior,  one  may  wonder  whether  even 
better  performance  can  be  obtained  with  more  components.  Table  4 
shows  the  performance  of  the  spectrogram  correlator  at  0 dB  SCR 
and  0 dB  SNR  for  five,  seven,  and  nine  principal  components.  The 
table  indicates  that  very  little  improvement  will  be  obtained  by 
using  more  than  seven  components  for  spectrogram  analysis. 


Table  4.  Performance  of  the  spectrogram  correlator  at 
0 dB  SNR  or  0 dB  SCR  for  5,  7,  and  9 principal 
components 


Classification  problem 

Number  of  components 

5 

7 

9 

Pr  [Misclassification, 

2.  5 x 7"  Al.  cyl. 
vs.  other  targets] 

in  noise 
(0  dB  SNR) 

0.  17 

0.  10 

0.093 

in  reverb. 

(0  dB  SCR) 

0.  32 

0.  27 

0.  25 

Pr  [Misclassification, 
cylindrical 
vs.  non-cylindrical 
targets] 

in  noise 
(0  dB  SNR) 

0.073 

0.051 

0.  048 

in  reverb. 

(0  dB  SCR) 

0.  14 

0.  14 

0.  13 

Pr  [Misclassification, 
metal 

vs.  non-metal 
targets] 

in  noise 
(0  dB  SNR) 

0.058 

0.  055 

0.052 

in  reverb. 

(0  dB  SCR) 

0.  12 

0.  12 

0.  12 
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1.6  Conclusion 

The  best  software  algorithm  to  date  is  a spectrogram 
correlator  that  uses  principal  component  analysis.  For  a single 
pulse,  this  processor  gives  a 0.  1 probability  of  misclassification 
for  a worst-case  signal-to-clutter  ratio  of  6 dB.  If  multiple  pulses 
are  used,  the  noise  component  of  the  measured  clutter  process  will 
be  reduced.  The  effect  of  the  clutter  itself  can  be  reduced  if  a 
target  can  be  viewed  from  a sequence  of  different  aspect. angles, 
since  different  bottom  reflections  will  be  encountered  at  each 
aspect,  and  the  clutter  echoes  will  be  more  noise-like  or  random. 

The  results  of  Figure  7a  indicate  that  a small  aluminum  cylinder  can 
be  distinguished  from  other  targets  even  when  the  aspect  of  the  cylinder 
is  varied.  Additional  pulses  and  "look"  angles  can  therefore  be  expected 
to  further  improve  performance  in  an  operational  environment. 

The  use  of  a spectrogram  correlation  process  for  detection 
of  range-distributed  targets  can  be  defended  from  three  different 
viewpoints: 

1.  It  works  better  than  anything  we  have  tried  thus  far. 

2.  A spectrogram  representation  is  apparently  a good 
approximation  to  the  representation  of  auditory  data 
that  is  passed  to  the  central  nervous  system  in  mammals 
(Volume  2,  Chapters  2.  2 and  2.  3). 

3.  Spectrogram  correlation  is  equivalent  to  optimum  Bayes 
or  Neyman- Pearson  detection  for  range-distributed 
targets  at  low  signal-to-noise  ratio  (Volume  2,  Chapter  2.4). 

The  spectrogram  correlator  as  an  optimum  detector  and  as  a 
model  for  auditory  signal  processing  is  discussed  in  Volume  2. 
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Appendices  for  Volume  1. 

APPENDIX  A:  FORTRAN  PROGRAM  FOR  CONSTRUCTION  OF  A 
SHORT- DURATION  SPECTROGRAM 

C 

C#****  MAIN  A L TCP  . <5 0 R A M .1 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c*#*** 

c 

c 

c 

c 

c***** 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SI IF  NO SC » SAN  DIEGO 

MACHINE  IJN I VAC  1.1. 10 

COMPILER  - FTN< ASCII ) FORTRAN  15  AUG  1978 

EXT ♦ libraries: 

ALTES * 

N*FTNL  IB. 

function: 


COMPUTE  SPECTROGRAMS  OF  TARGET  ECHO  PLUS  NOISE  SAMPLES 


PROCESS  CONTROL  PARAMETERS: 


PARAMETER 

NTARG 

MTARGd  ) 


SCI 


NSMP 


USAGE 

NUMBER  OF  TARGET  ECHOS  TO  BE  PROCESSED 
<D=38»  R-l-38) 

RECORD  POSITION  IN  RAN  ACCESS  FILE  OF 
I'TH  TARGET  <D=I»  R=l-38) 

S I GNAL  COMPOS I T I ON  INDEX 
I F < SC I - 1 ) signal:  SF*ECHO+NOISE 
IF<SCI=0)  SIGNAL.;  ECHO 

IF ( SCI= -1 ) signal:  noise 

NUMBER  OF  NOISE  SAMPLES  TO  BE  USED 
PER  TARGET 


RNIO  LABEL  OF  FIRST  NOISE  SAMPLE.  SEQUENTIAL 

ACCESS  IS  USED  THEREAFTER.  (D--1  ) 


NSNR 


NUMBER  OF  SIGNAL  TO  NOISE  RATIOS  TO  USE 
PER  TARGET  PER  NOISE  SAMPLE  <D~6>  R=1 - 6 ) 


AS NR ( I ) 


VALUE  OF  THE  I'TH  SIGNAL  TO  NOISE  RATIO 
( D “ S E E DATA  ST  AT  E MEND 
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APPENDIX  A (continued) 


C LDSPLA 

C 

c 

LOGICAL  VARIABLE  IND  I CATES  WHETHER  SPEC  I RO 
GRAM  IS  TO  BE  DISPLAYED  <=T)  OR  NOT  ( =F ) 

( D=F ) 

t 

c*****  FILES 

C FILE 

USE 

C TEF/V 

C NF/3 

C FMF/ 8 

C RF/10 

C CARDS/5 

C PRINTER/6 

- INPUT  TARGET  ECHO  FILE 

- INPUT  NOISE  FILE 

- INPUT  MATCHED  FILTER  FILE 

OUTPUT  RESULTS  FILE 

- INPUT  CARD  FILE 

- OUTPUT  PRINT  FILE 

C 

C 

C*****  DATA  STRUCTURE  ELEMENTS 
C 

c 

DIMENSION  ASNR  ( 1 0 ) . ECHO  < 1 024  ) , EPN  < 1 024  ) . EN(3  ( 200 ) . W ( 32 ) , 

1 SF'CG  < 28 . 1 6 ) , I SPCG  (28.16) 

I NTEGER  MTARG ( 38 ) * SC I . TEE  * FMF » RF 
REAL  NOIS< 1024 ) 

COMPLEX  SMF (512) » CEPN ( 1024 ) . F ( 32 ) 

LOGICAL  LDSPLA 
C 

C*****  NAMELIST . DATA  STATEMENTS » FORMATS » ETC. 

C 

C 

NAMEL IST/CARDS/NTARG » MTARG .SCI » NSMP » RN 1 0 . NSNR . ASNR . LDSPLA 
DATA  ( ASNR ( I ) » 1=1 . 6 ) /100 . . 20 . . 10 . . 4 . . 2 . . 1 . / 

1 FORMAT <///5X. 'ERROR  - TARGET  NUMBER  '.13.'  REQUESTS  ' , 

*' RECORD  '.I3> 

2 FORMAT (//10X. 'RUN  SPECIFICATION'// 

*5X . ' NO . OF  ECHOS ' . T40 . I 1 0 . T55 . ' NTARG ' / 

*5X.  'SIGNAL.  COMPOSITION  INDEX'  .T40. I10.T55.  'SCI '/ 

*5X. 'INITIAL  NOISE  SEQUENCE  TAG ' . T40. F.10 . 0 > T55 . 'RNIO'/ 

*5X » ' NO . OF  SAMPLES  PER  ECHO ' . T40 . 1 1 0 . T55 . ' NSMP ' / 

*5X . ' NO  . OF  SNRS  PER  SAMPLE  ' . T40 . 1.10 » T55 . ' NSNR  ' ) 

C 

1.1  - SET  PROCESS  CONTROL  PARAMETER  X DEFAULT  VALUES f 

C READ  CONTROL  PARAMETERS » READ  SIGNAL  MATCHED 

C FILTER 

C 
C 

NTARG “38 
DO  100  1=1.30 
MTARG ( I ) = I 
1.00  CONTINUE 
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APPENDIX  A (continued) 

3Ct  .1 
NF::-3 

El-  1 0 

TEF  :9 

FME-8 
RNIO-1  , 

N SNR-6 

LDSPLA-.  FALSI.  . 

RNI--0 » 

READ ( 5 • CARDS > END  I 1 0 ) 

1 10  CONTINUE 

READ  ( Ri I ) < SHF  ( I ) * ]>1 . 51.2  ) 

IF  (SCI  . EQ.  1.  'GO  It)  1.40 
IF  (SC  I.  EG.  1)00  TO  .120 
A SNR  ( 1 ) ■■■■■  I .1  '.() 

NSMP=1 
GO  TO  1.30 
120  CONTINUE 
AS  NR  ( 1 > -0 . 

NTARG-"  1 
130  CONTINUE 
NSNR--1 

140  CONTINUE 

UR I TE ( 6 t 2 ) NTARG r SC I , RN I 0 , NSMP , NS NR 

C 

C#****  1.2  BEGIN  TARGET  ECHO  I OOP 

C 

C INITIAL  CONDITION  - ITARG  -1 

C TERMINATION  CONDITION  - ITARONTARG 

C RETURN  POINT  LABEL  - 200 

C EXIT  POINT  LABEL  - 620 

C 

c 

DEE  I NE  f I LE  9 ( 38, 1 024  > IJ . N? ) 

ITARG=1 
200  CONTINUE 

IF  < I T ARG • GT . NTARG ) GO  TO  620 
C 

C#****  2.1  READ  TARGET  ECHO  IE  REQUIRED 

C 

C 

IF ( SCI .GE. 0)00  TO  210 
NER-0 

DO  205  II .1024 
ECHO  ( I ) ::::0 . 

205  CONTINUE 
GO  TO  220 
210  CONTINUE 


[ 


i 
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APPENDIX  A (continued) 

NER-TM  T ARG  ( I T ARG  ) 

1 r ( NER . LE . 38 . AND . Nt  R . GE  . I ) GO  TO  2 1 5 

WRITE (A, 1 ) 1 TARG » NE  R 

STOP 

215  CONTINUE 

READ (TEF' NER) ECHO 
220  CONTINUE 
C 

(;.*****  2.2  BEGIN  NOISE  SAMPLE  LOOP 

C 

C INITIAL  CONDITION  - ISMP=1 

C TERMINATION  CONDI IION  - ISMP-NSMP 

C RETURN  POINT  LABEL  - 300 

C EXIT  POINT  LABEL  - 610 

C 
C 

ISMP=1 

300  CONTINUE 

IF< IShP.GT .NSMP )G0  TO  610 
C 

<;♦****  3.1  READ  NOISE  SAMPLE  IF  REQUIRED 

C 

C 

IF (SCI ,EQ.O)GO  TO  310 

CALL  RDNOIS ( NF . NOIS » RNI r RNIO ) 

GO  TO  330 
310  CONTINUE 

DO  320  1=1 r 1024 
NOISd  )=0. 

320  CONTINUE 
330  CONTINUE 
C 

C***#*  3.2  BEGIN  SIGNAL  TO  NOISE  RATIO  LOOP 

C 

C INITIAL  CONDITION  - ISNR=1 

C TERMINATION  CONDITION  - ISNR>NSNR 

C RETURN  POINT  LABEL  - 400 

C EXIT  POINT  LABEL  - 600 

C 
C 

ISNR-1 

400  CONTINUE 

IF ( I SNR . GT . NSNR ) GO  TO  600 
SNR=ASNR( ISNR) 

C 

4.1  MIX  ECHO  AND  NOISE  BY  SNR  AND  PASS  THRU 
C FILTER  MATCHED  TO  SIGNAL 

C 
C 
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APPENDIX  A (continued) 


I Nh'SNI 

II  (SCI  ,1  U o < INK  I * 

I A I I I NM  I * < L < HU  . N(  1 1 S , JR  , EPN  » 

MO  410  1-1*1024 

CF. I'M ( 1 1 =CMPLX ( EPN < I ) * 0 . > 

410  CONI INUE 

CAL  I RF  F I ( CEF'N  * 1 024  *10*0) 
nil  420  1 = 1*512 

C-F  I 'M  < r )=cn  'N(  l >*::;mf  < i > 

420  CONTINUE 
C 

C****#  4.2  LOW  PASS  FILTER  RESPONSE  TO  5 A KHZ  AND 

C RETURN  TO  TIME  DOMAIN  AND  SAMPLE  RESPONSE 

C 

DO  430  1 = 11  A*  1.024 
CEF'N ( I >=  CMPI  X(().  *0.  > 

430  CONTINUE 

CEF'N  (114  ) =CEPN  (114)  KCMF'LX  ( . 5886  * 0 . ) 

CEF'N  (115)  = CEF'N  (115)  #CMF'LX  ( . 1 0A5 , 0 . ) 

CALL  REF  T (CEF'N*  1024*  10*  1 ) 

K=0 

DO  440  1=4*1024*4 
K=KF1 

EPN  ( K ) = REAL  ( CEF'N  ( I ) ) 

440  CONTINUE 
C 

C*****  4.3  SMOOTH  ENERGY  FUNCTION  AND  LOCATE  MAX 

C 

C 

Sl!M2=0 . 

DO  450  1=1 *1A 
SUM2=SUM2+EPN ( I ) ** 2 
450  CONTINUE 

ENG  ( 1 ) -=SUM2 
EMAX=ENG( 1 ) 

KMAX=t 

DO  470  1=2*256 
K=I f 15 

IF ( K . GT . 25A > K=K-25A 

S U M 2 = S U M 2 + E P N < K > * * 2 •-  EPN  (I- I )*  * 2 

I F ( SUM2 . t.E . EM  AX ) GO  TO  4A0 

EMAX=SUM2 

KMAX=I 

4 AO  CONTINUE 

ENG ( I ) =SUM2 
470  CONTINUE 
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APPENDIX  A (continued) 


C*#***  4.4  CONSTRUCT  SPECTROGRAM  J WINDOW  DATA  IN 

C TIME  domain;  transform  overlapped  segments 

c 

c 

DO  500  1=1,15 
X = I -0.5 

U ( I ) = . 5*  < 1 . +COS (3.141 S9*X/7 . 5 ) ) 

500  CONTINUE 

DO  510  1=16,32 
ENG  < I ) =0 . 

510  CONTINUE 

DO  550  IND=1 » 28 
IS=KNAX+  < IND-l ) #3-32 
DO  515  1=1,15 
J=IS+I 

IF< J.LT.l) J=J+256 
IF < J.GT.256) J=J-256 
ENG  < I > =U ( I ) *EPN ( J ) 

515  CONTINUE 

DO  520  1=1,32 
F ( I > =CMPLX < ENG ( I ) , 0 . ) 

520  CONTINUE 

CALL  RFFT<F ,32,5,0) 

DO  550  1=1,16 

SPCG  < IND, I >=F< I > *CON JG (F ( I ) ) 

550  CONTINUE 
SMAX=0. 

DO  560  N=1 » 28 
DO  560  1=1,16 
SMAX=AMAX1 ( SPCG ( N , I ) , SMAX ) 

560  CONTINUE 

DO  570  N=1 , 28 
DO  570  1=1,16 

SPCG (N, I ) =SPCG < N , I ) *99 . /SMAX 
I SPCG ( N , I ) =SPCG  < N , I ) 

570  CONTINUE 
C 

C****»  4.5  WRITE  RESULTS  TO  OUTPUT  FILE  AND  DISPLAY 

IF  LDSPLA-T 


WRITE  (RF  )NER,SNR,RNI , ( ( SF’CG(  J,  I ) , J=l,28,2)  , 1 = 1,16) 
IF  ( L.DSPLA ) CALL  DSPLA  ( NER , SNR , RNI , I SPCG ) 

ISNR=ISNRH 
GO  TO  400 
C 

(;;*****  3.3  END  SIGNAL  TO  NOISE  RATIO  LOOP 
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APPENDIX  a (concluded) 


600  CONTINUE 

ISMP-IShP+l 
GO  TO  300 
C 

C*****  2.3  ENI)  NOISE  SAMPLE  LOOP 

C 

c 

610  continue 

I TARG= I TARG+ 1 
GO  TO  200 
C 

C*****  1.3  END  TARGET  ECHO  LOOP 

C 

c 

620  CONTINUE 
ENDFILE  RE 
ENDF1I  E RF 
END 


s 
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APPENDIX  B:  FORTRAN  PROGRAM  FOR  GENERATING 
FM-GRAM 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c**### 

c 

c 

c 

c 

c 

c##*## 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


21  AUG  1978 


MAIN  AI.TES . FGRAM1 


SITE:.  NO  SC . SAN  m EGO 

MACHINE  UNI VAC  1110 

COMPILER  FTN(ASCII)  FORTRAN 
EXT.  LIBRARIES t 

ALTES. 

N*FTNLIB. 


function: 

FORMS  A MODIFIED  SPECTROGRAM  THAT  IS  BASED 
UPON  MATCHED  FILTERS  WITH  DIFFERENT  LOG 
CHIRP  FACTORS 

PROCESS  CONTROL  PARAMETERS 
PARAMETER  USAGE 


NTARG 


MTARG ( I ) 


SCI 


NSMP 


RNIO 


NSNR 


ASNR ( I ) 


.DSPLA 


NUMBER  OF  TARGET  ECHOS  TO  BE  PROCESSED 
( D=38 f R-l-38) 

RECORD  POSITION  IN  RAN  ACCESS  FILE  OF 
I'TH  TARGET  (D~I»  R=l~38) 

SIGNAL  COMPOSITION  INDEX 
IF" ( SC  1 = 1 ) signal:  SF#ECHOf noise 
IF (SC 1=0)  signal:  echo 
if<sci=-i>  signal:  noise 

NUMBER  OF  NOISE  SAMPLES  TO  BE  USED 
PER  TARGET 

LABEL  OF  FIRST  NOISE  SAMPLE.  SEQUENTIAL 
ACCESS  IS  USED  THEREAFTER.  <D=1> 

NUMBER  OF  SIGNAL  TO  NOISE  RATIOS  TO  USE 
PER  TARGET  PER  NOISE  SAMPLE  ( D=A • R~l-A) 

VALUE  OF  THE  I'TH  SIGNAL  TO  NOISE  RATIO 
(D-SEE  DATA  STATEMENT) 

LOGICAL  VARIABLE  INDICATES  WHETHER  SPECTRO- 
GRAM IS  TO  BE  DISPLAYED  <=T>  OR  NOT  <=F) 

Oi  l ) 
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APPENDIX  B (continued) 


' 


i ; 


C****«  FILES 

C 

C 

FILE 

USE 

l, 

C 

TEF/9 

- INPUT  TARGET  ECHO  FILE 

f 

NF/3 

- INPUT  NOISE  FILE 

c 

RF/10 

OUTPUT  RESULTS  FILE 

c 

CARDS/5 

- INPUT  CARD  FILE 

c 

c 

PR INTER/ A 

- OUTPUT  PRINT  FILE 

C 

C*****  DAT  A STRUCTURE  ELEMENTS 

C 

DIMENSION  ASNR  < 1 5 ) * ECHO ( 1024) , RESP (7,256) *EPN( 1024) , IRESP ( 7 1 100 ) , 


1RESP1 (7, 100 ) 

INTEGER  MTARG ( 38 > , RE , SCI y TEE 
REAL  NO IS (1024) 

COMPLEX  U ( 7 y 1 20 ) * V < 256 ) , DD ( 256 ) 

LOGICAL  LDSPLA 
C 

C*****  NAMELIST,  DATA  STATEMENTS,  FORMATS,  ETC. 

C 

c 

DEFINE  FILE  9(38, 1024, U,N9> 

NAMELIST  /CARDS/NTARG , MTARG , SC I , NSMP , RNIO , NSNR , ASNR , LDSPLA 
DATA  ( ASNR ( I ) , 1-1 , 6 )/100 . , 20 . , 10 . , 4 . , 2 . , 1 . / 

1 FORMAT <///5X» 'ERROR  - RECORT  NUMBER  ',12,'  REQUESTS  ', 

* ' RECORD  ',13) 

2 FORMAT <//10X, 'RUN  SPECIFICATION'// 

*5X,'N0.  OF  ECHOS'  ,140,  II  0,T55,  'NT  ARGV 

*5X,  'SIGNAL..  COMPOSITION  INDEX ' , T40 , 1 10  , T55 , ' SCI ' / 

*5X , ' IN I T I AL  NO I SE  SEQUENCE  TAG ' , T 40 , F 10.0, T55 , ' RN 1 0 ' / 
*5X,'N0.  OF  SAMPLES  PER  ECHO ' , T40 , I 10 , T55 , ' NSMP ' / 

*5X , ' NO . OF  SNRS  PER  SAMPLE ' , T40 , I 1 0 , T55 , ' NSNR ' / 

#5X r 'S/N  RATIOS' ,T40,E10.4,T55, 'ASNR(I)  '/ 

*14(140, E10.4/)  ) 

C 

C*****  SET  PROCESS  CONTROL  PARAMETERS  AND  DEFAULT  VALUES 
C 

C READ  CONTROL  PARAMETERS 

C 

NTARG-38 
DO  100  1=1,38 
MTARG (I >=I 
100  CONTINUE 
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APPENDIX  B (continued) 


I 


sc.;  i=i 

TEF  ~9 
NF=3 
RF-10 
RN 1 0*  :l . 

NS NR « 6 

LDSPLA-- . I:  AI...SE . 

EN I ”0 . 

READ ( 5 . CARDS » END” 1 1 0 ) 
l 1.0  CONTINUE 

IF (SCI .EG. 1 )G0  TO  140 
IF(SCI.EG,  1)00  TO  120 
A SNR ( 1 ) ” 1 . E30 
NSMP- 1 
GO  TO  130 
120  CONTINUE 
ASNR ( 1 ) ”0 . 

NTARG-l 
130  CONTINUE 
NSNR~ 1 

140  CONTINUE 

WR I TE (6*2) NTARG * SC  I . RN 1 0 . NSMP » NSNR , ( ASNR (J), J- 1 . NS NR ) 

C 

C*****  COMPUTE  ' MATCHED ' FILTERS  WITH  RANGE  OF  CHIRP  FACTORS. 

C 

C 

S-47 . 

FK:=1  • 2 
FC=- . 25 

WN-l . / ( 2 . *ALOG ( FK ) ) 

W=0 . 

DO  160  N"1 . 7 
DO  160  1=1. 128 

IFd.LT.10  .OR.  I. GT. 127)00  TO  150 
X2~  < 1-1 . >/S 
Yl-ALOG < X2 ) 

Y5:=;  1 , / X2 

IFd.EQ.10  .OR.  I ,EQ.  127) Y5~0« 1065*Y5 
I F < I . EQ . 1 1 . OR . I . EO . 1 26 ) Y5”0 , 5886* Y5 
CHRP=N-3 , 

L~ 2 . *WN*6 , 2832* Y 1 *FC*CHRP 
CN-COS ( 2 ) 

SN=SIN<Z> 

Y6=Y5*CN 

Y7=Y5*SN 

U ( N . I ) kCMPLX < Y6  * ■ Y7 ) 

GO  TO  160 
150  CONTINUE 
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APPENDIX  B (continued) 


U t N » 1 ) -CMPLXvO. .0. ) 

1 AO  CONTINUE 
C 

C*#***  BEGIN  TARGET  ECHO  LOOP 

C 

C INITIAL  CONDITION  •••  ITARG=1 
C TERMINATION  CONDITION  - ITARG>NTARG 
C RETURN  POINT  LABEL  - 200 
C EXIT  POINT  LABEL  - A20 
C 
C 

ITARG-1 
200  CONTINUE 

IF < ITARG . GT . NTARG ) GO  TO  A20 
C 

C***»*  READ  TARGET  ECHO  IF  REQUIRED 

C 

C 

IF (SCI .GE*0)G0  TO  210 
NER=0 

DO  205  1=1,1024 
ECHO ( I ) =0 . 

205  CONTINUE 
GO  TO  220 
210  CONTINUE 

NER=MTARG( ITARG) 

IF (NER.LE . 38 .AND .NER .GE . 1 )GQ  TO  215 

WRITE  < 6 » 1 ) ITARG  , NER 

STOP 

215  CONTINUE 

READ (TEF' NER) ECHO 
220  CONTINUE 
C 

C*#***  BEGIN  NOISE  SAMPLE  LOOP 
C 

C INITIAL  CONDITION  - ISMP=1 
C TERMINATION  CONDITION  ~ 1SMP>NSMP 
C RETURN  POINT  LABEL  - 300 
C EXIT  POINT  LABEL  - 610 


ISMP=1 

300  CONTINUE 

IF  < ISMP . GT  » NSMP ) GO  TO  610 

C 

C*****  READ  NOISE  TARGET  IF  REQUIRED 

C 

C 
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APPENDIX  B (continued) 


IF<SCI.EO.O)GO  10  310 

CALL  RDNO I S < NF » NO IS, RN I ,RNIO> 

GO  TO  330 
310  CONTINUE 

DO  320  1=1,1024 
NOIS< I )=0. 

320  CONTINUE 
330  CONTINUE 
C 

C***#*  BEGIN  SIGNAL  TO  NOISE  RATIO  LOOP 
C 

C INITIAL  CONDITION  - ISNR=1 
C TERMINATION  CONDITION  - ISNR>NSNR 
C RETURN  POINT  LABEL  - 400 
C EXIT  POINT  LABEL  - 600 
C 

ISNR=1 

400  CONTINUE 

IF<  ISNR'.GT  .NSNR)GO  TO  600 
SNR=ASNR ( ISNR ) 

C 

C*****  MIX  ECHO  AND  NOISE  BY  SNR 
C 

c 

TNR=SNR 

IF<SCI.EG.0)TNR=1. 

CALL  ENMIX ( ECHO , NQIS  > TNR , EPN ) 

C 

C#****  COMPUTE  TRANSFORM  OF  SAMPLED  SIGNAL 

C 

C 

DO  410  1=1,256 
K=4*I 

DD ( I ) =CMPLX ( EPN ( K ) , 0 • ) 

410  CONTINUE 

DO  420  1=1,256 
EPN ( I ) =REAL ( DD ( I ) ) 

420  CONTINUE 

CALL  RFFT  < DEi , 256 ,8,0) 

C 

C***#*  COMPUTE  'MATCHED ' FILTER  RESPONSES. 
C 

c 

DO  450  N=1 , 7 
DO  430  1=1,128 
V(I)=U(N,I)*DD(I) 

V< 1+128 )=CMPLX(0. ,0. ) 

430  CONTINUE 
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APPENDIX  B (continued) 


CALL  REFT  < V » 256 » 8 » t ) 

DO  440  1=1,256 
RESP < N • I ) =CADS ( V ( I ) )**A 
440  CONTINUE 
450  CONTINUE 
C 

C***#*  NORMALIZE  RESPONSES 
C 
C 

RMAX=0 . 

IMAX=0 

DO  460  N=1 » 7 
DO  460  1=1,256 
IF < RESP < N , I ) . LE . RMAX ) GO  TO  460 
IMAX=I 

RMAX=RESP  < N » I ) 

460  CONTINUE 

DO  470  N=1 » 7 
DO  470  1=1,256 

RESP (Nfl) =RESP (N , I ) *99 . /RMAX 
470  CONTINUE 

IBEG=IMAX~49 

IEND=IMAX+50 

K=1 

DO  480  I = IBEG , IEND 
DO  475  N=1 , 7 
J=I-8*(N~4) 

IF(J.LT.l) J= J+256 
IF(J.GT . 256 ) J= J-256 
RESP1 <N»K)=RESP(N» J) 

I RESP  < N f K ) =RESP 1 ( N » K ) 

475  CONTINUE 
K=K+ 1 

480  CONTINUE 

C 

C*****  WRITE  RESULTS  TO  OUTPUT  FILE  AND 
C DISPLAY  IF*l.DSPLA=T 

C 

WRITE < RF ) NER t SNR t RNI » ( ( RESP1 ( N » I ),  1 = 1 1 100»2) > N=1 f 7) 
IF ( L DSP LA ) CALL  DSPLA ( NER » IRESP ) 

!SNR=ISNR+1 
GO  TO  400 
C 

c******  END  SIGNAL  TO  NOISE  RATIO  LOOP 

C 

C 

600  CONTINUE 

ISMP=ISMP+1 
GO  TO  300 
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APPENDIX  B (concluded) 


C 

(.;**#**  END  NOISE  SAMPLE  LOOP 
C 

c 

A 10  CONTINUE 

I TAEG~ I TARGE 1 
GO  TO  200 
C 

C*****  END  TARGET  ECHO  LOOP 
C 

c 

620  CONTINUE 
END 


